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Abstract 

A magnetohydrodynamic model is constructed for a cylindrical jet immersed in an ex- 
ternal uniform magnetic field. It is shown that, as in the force-free case, the total electric 
current within the jet can be zero. The particle energetics and the magnetic field structure 
are determined in a self-consistent way; all jet parameters depend on the physical conditions 
in the external medium. In particular, we show that a region with subsonic flow can exist in 
the central jet regions. In actual relativistic jets, most of the energy is transferred by the elec- 
tromagnetic field only when the magnetization parameter is sufficiently large, a > 10 6 . We 



o 



also show that, in general, the well-known solution with a central core, B z = B /(l+m 2 /w 2 c 
cannot be realized in the presence of an external medium. 

o 



1 Introduction 



o 

of compact astrophysical objects. Indeed, jets are observed in most compact sources, ranging 



X 



The formation mechanism of jets is a key issue in the study of the magnetospheric structure 



from active galactic nuclei (AGNs), quasars, and radio galaxies |j| to accreting neutron stars, 
solar-mass black holes (SS 433, X-ray novae) 0, and young stellar objects 0. Moreover, jets 
have also been recently discovered in young radio pulsars [|], ||. At the same time, in most 
studies devoted to the magnetohydrodynamic (MHD) model of such objects || |7|, [8], [II], [11 



m 



which the formation of jets is coupled with the attraction of longitudinal currents flowing in the 
magnetosphere, the attention was focused on intrinsic collimation in the sense that the effect of 
the external medium was assumed to be marginal. However, such a situation is possible only for a 
nonzero total current I flowing within the jet |T2|], so the question of its closure in the outer parts 



of the magnetosphere arises. On the other hand, the longitudinal current is often constrained 
by the regularity condition at the fast magnetosonic surface, which by no means always leads to 



the sufficiently large longitudinal currents required for collimation | 13fl . In other words, with the 



exception of the force-free case fll4| , as yet no working model of a jet in which, on the one hand, 



the total electric current would be zero and, on the other hand, the total magnetic flux in 
the jet would be finite, has been constructed. However, the force-free approximation, in which, 
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by definition, the particle energy density is disregarded, does not allow the fraction of energy 
transferred by the outflowing plasma to be determined. 

At the same time, the question of collimation cannot be solved in isolation from the external 



conditions (see, e.g., |L5], |16|). In particular, this is clear even from a popular example of the 
magnetosphere of a compact object with a monopole magnetic field, because for any arbitrarily 
weak external regular magnetic field, the monopole solution (for which the magnetic field falls off 
as r~ 2 ) cannot be extended to infinity. Moreover, as is well known from an example of moving 
cosmic bodies, such as Jupiter's moons [|T7|] or artificial Earth satellites [1^, [HJ, as well as radio 



pulsars [pOR , the external magnetic field can serve as an effective transfer link, which occasionally 
determines the general energy losses of the system. For this reason, constructing a consistent 
magnetospheric model for compact objects immersed in an external magnetic field is, in our view, 
of undeniable interest, especially since, as was noted above, such a jet model was previously 



constructed in the force- free approximation [14 



Undoubtedly, the existence of an external regular magnetic field in the vicinity of compact 
objects is largely open to question. The regular magnetic field in our Galaxy, i.e., the field that is 
constant on scales comparable to the sizes of our Galaxy, is known to be 

5 ext ~10- 6 G, (1) 

and essentially matches the random magnetic-field component, which varies even on scales of 
several parsecs [21] . However, if the collimation is assumed to be actually produced by an external 
magnetic field, it becomes possible to estimate the jet radius. Indeed, assuming the magnetic field 
in the jet to be similar to the external magnetic field (|l|), we obtain from the condition for the 
conservation of magnetic flux 

/ B- \ 1/2 

where R and B[ n are the radius and magnetic field of the compact object, respectively. For 
example, for AGNs (B in ~ 10 4 G, R ~ 10 13 cm), we have 

rj ~ 1 pc, (3) 

which corresponds to the observed jet radii One might expect such a picture to be also 
preserved for an external medium with pressure P ~ B 2 /8ir; therefore, it seems of interest to 
consider the internal structure of a one-dimensional jet immersed in an external uniform magnetic 
field. However, a discussion of the more realistic case of a medium with pressure is beyond the 
scope of this study. Nor do we discuss the collimation itself but only consider the internal structure 
of observed one-dimensional jets. This issue has become particularly urgent because of the new 
possibilities offered by space radio interferometry, which enables the internal structure of such 
jets to be resolved. The effect of an external medium on the internal structure of relativistic jets 
in the MHD model discussed here was previously studied only by Appl and Camenzind Jl5|, [L6[] . 
They considered only a special case with a constant angular velocity of the plasma, in which the 
solution with a zero total electric field flowing inside the jet could not be constructed. As we 
show below, it is for the case of an angular velocity decreasing toward the jet periphery (which, 
incidentally, is typical of all models with a magnetic field passing through the accretion disk) that 
the solution with finite magnetic flux \l/o an d zero total current /(^o) — can be constructed. 
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On the other hand, many authors |22|, [23|, |24|, obtained a universal solution with a central 
core for a cylindrical jet: 

B * = 14.5/^ ' (4) 

where, in the relativistic case, 

w e =- (5) 

is the size of the central core, Q is the angular velocity of the compact object, and 7 is the 
characteristic Lorentz factor of the outflowing plasma. As can be easily seen, such a solution results 
in a rapid falloff of the poloidal field B z oc vo~ 2 far from the rotation axis w ^> vj c . However, 
this solution is in conflict with the force-free approximation, in which the poloidal magnetic field 



remains essentially constant [lj]. Indeed, when the energy density of the electromagnetic field 
exceeds appreciably the plasma energy density (and it is this case that was considered), it would 
be natural to assume that the internal jet structure must be similar to the force-free one. 

The examples given above show that a more detailed study with allowance for all possible 
solutions is required even in the simplest case of a one-dimensional cylindrical jet considered in 
terms of ideal magnetohydrodynamics. Our study aims at a consistent investigation of this issue. 

Several features that we use when studying the structure of relativistic jets typical of AGNs 
and radio pulsars should be immediately noted. First of all, the jet radius r } in all real cases proves 
to be considerably larger than the light-cylinder radius = c/Q. This implies that, when the 
internal structure of jets is investigated, the corresponding equations must be written in complete 
relativistic form. On the other hand, the gravitational forces can be disregarded in them far from 
the compact object. Finally, for simplicity, we consider below a cold plasma, which is justifiable 
because the thermal processes in the magnetospheres of radio pulsars play no crucial role. As 
for the jets from AGNs, this approximation is applicable here in those magnetospheric regions in 
which the plasma density is low. In any case, this is true for the field lines passing through the 
surface of a black hole. 

The one-dimensional solutions describing collimated jets are obtained in Sect. 2; analytic and 
numerical solutions for the basic physical quantities characterizing the structure and physics of jets 
are given in Sect. 3. The problem is solved in straightforward statement; i.e., all jets characteristics 
are determined by a set of parameters in the compact source and, most importantly, by the physical 
conditions in the external medium. As a result, we have found the conditions under which most of 
the energy in actual relativistic jets must be transferred by the electromagnetic field, while a region 
with subsonic flow exists in the central jet regions. We also show that the solution with a central 
core (|4]) and @ cannot be realized in an external magnetic field. Finally, some astrophysical 
implications of the theory developed for one-dimensional jets are discussed in Sect. 4. 



2 Basic equations 

Let us consider the structure of a one- dimensional jet where all quantities depend only on radius 
w; in what follows, the temperature of the matter is assumed to be zero, and c = 1. As in the 
general axisymmetric case, it is convenient to describe the magnetic-field structure in terms of 
magnetic-flux function ^f(w), which is related to the longitudinal magnetic field by 

Zttzu aw 
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Accordingly, it is convenient to write the toroidal magnetic field, the electric field, and the 
4-velocity vector of the matter as 

B v (w) = (7) 



2tt aw 

u = — B + ■yfl-pwe^, (9) 
n 

where I{wq) is the total current within w < w . In the case of a cold plasma, at the cylindrical 



magnetic surfaces ^ = const, four "integrals of motion" |26| can be introduced, which should be 
considered precisely as functions of magnetic flux ^ in the most general statement. These are 
primarily fip^) and ry(^/) in the definitions (|8|) and (Q), as well as the z-component of angular 
momentum L(^f) = I/2tt + fir/zuu^ and the energy flux E(ty) = Q F I/2n + "f^rj. Here, /x is the 
relativistic specific enthalpy, which is equal to the mass of particles for a cold plasma. The specific 
form of the integrals of motion must be determined from boundary conditions in the compact 
source and from critical conditions at the singular surfaces. 

As a result, the equilibrium equation for magnetic surfaces far from gravitating bodies (Grad- 
Shafranov's equation) can be written as (see, e.g., f27j) 




^«V« F (v*)^ + ^A^_^*§!). , (io) 



wdw) v ' d$> w 2 2M 2 d^\AJ M 2 d^ 
where 

G = zu 2 e 2 + M 2 L 2 -M 2 zu 2 E 2 , 
A = 1 - Q 2 F w 2 - M 2 , 

e(#) = - ft F (*)£(*)- 
Here, M 2 = u 2 /u 2 is the square of the Mach number with respect to the Alfven velocity u a = 
5 2 /(47m/i) 1 / 2 , and the derivative d/d^f acts only on the integrals of motion. The remaining jet 



parameters are given by the well-known algebraic relations (see, e.g., [27]) 

I L- VL F w 2 E 



2ir ~ 1 - tt 2 F w 2 - M 2 ' 



1 (E-n F L)-M 2 E 
7 "~ m l-tt 2 ¥ w 2 -M 2 ' { ' 

1 (E - n F L)n F w 2 - LM 2 

Uip ~ w~^q 1 - Q 2 ™ 2 - M 2 " ( } 

Equation ( [TDD contains four integrals of motion; this equation has no singularity at the fast 
magnetosonic surface, because it depends only on coordinate w. As for the Alfven surface, A = 0, 
the problem of the boundary conditions generally requires a further study beyond the theory of 
ideal magnetohydrodynamics. At the same time, for the fairly large currents I ~ Iqj considered 
here, a solution continuous across this surface can always be constructed by a small change in the 
integrals of motion near the Alfven surface, $ ~ \&a- Consequently, equation (|I0D requires six 
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boundary conditions. These boundary conditions primarily include the external uniform magnetic 
field 

B z { ri ) = B ext , (14) 
and the regularity condition at the magnetic axis w — > 

m{w) -> Cw 2 . (15) 

In addition, all four integrals Qp, E, L, and rj must be specified. As for the remaining quantities 
characterizing the flow, such as the jet radius rj and the outflowing plasma energy, they must be 
determined as a solution of the problem formulated above. Similarly, the solution of the problem 
must also give an answer to the question of whether the flow in the jet is supersonic. 

Let us now consider the determination of the integrals of motion in more detail. It would be 
natural to assume that, at the jet boundary where there is no longitudinal motion of the matter, 
all four integrals of motion become zero 

fi F (^o) = 0, £(* ) = 0, £(* )=0, 77(*o) = 0. (16) 

Here, \l/o is the finite total magnetic flux concentrated in the jet. This case corresponds to the 
absence of tangential discontinuities at the jet boundary; according to (|TTD, the total electric 
current within the jet is automatically equal to zero. 

We use the integrals of motion f2 F (\I>), L(\&) = I(^)/2tt, and E($) = f) F derived by 
Beskin et al. [B^j for the force-free magnetosphere of a black hole, which satisfy the conditions ( ]T5| ) 
and, consequently, can be directly used to study the jet structure. The only but very important 
change here is the fact that, for a finite magnetization parameter a J29j, 



, = (17) 
8tt 2 /^(O)' V ; 

which tends to infinity in the force-free approximation, the particle contribution must be added 
to the energy integral E(^f), because the energy flux of the electromagnetic field near the rotation 
axis must inevitably vanish. As a result, we have with accuracy up to ~ a -1 <C 1 



2J1 - 

n p (*) = v . n r (o), 

1 + Jl - ^/^ 



1 , 1 - ^/^ 

m = 7n i ^(°)^ ( 19 ) 



E(^)= linf ir ] + Q F ^)L^). (20) 
Below, we assume, for simplicity, that 

7i n = const. (21) 

We emphasize that 7; n in expression ( ^0|) has the meaning of the injection Lorentz factor in 
the region of the compact object and it is not equal to the Lorentz factor of the jet particles. 
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Thus, we see from the formula for the energy flux E(ty) that the contribution by the electro- 



magnetic field becomes dominant only at ^ > ^/j n , where 



a 



(22) 



At low values of most of the energy is transferred by the relativistic particles; as directly 
follows from relation fl2"D|), their Lorentz factor is constant and equal to their initial value 7^. As 
for the integral T](^f), the particle-to-magnetic flux ratio, we chose it in the form 



77(tf) = 770(1 - tf/#o), 



(23) 



which satisfies the condition (TO). 

We emphasize that the very possibility of using the integrals of motion obtained by analyzing 
the inner magnetospheric regions, is not trivial. Indeed, the flow outside the fast magnetosonic 
surface is completely determined by four boundary conditions at the surface of a rotating body. 
At the same time, a one-dimensional flow can be produced by the interaction with the external 
medium, which gives rise (see, e.g., |]3"0|) to perturbations or shock waves propagating from "acute 
angles" and other irregularities. Therefore, in regions where the conditions for the validity of 
ideal magnetohydrodynamics are violated, a significant redistribution of energy E and angular 
momentum L is possible (e.g. a part of them can be lost via radiation). Nevertheless, we assume 
here, for simplicity, that the integrals of motion Efif) and Lfif), functions of flux ty, remain 
exactly the same as those in the inner magnetospheric regions. 

In the one-dimensional case we consider, it is convenient to reduce the second-order equa- 
tion fllCf ) to a set of two first-order equations for ty(m) and M 2 {m). Multiplying equation ([H]) by 
2A(d^ /dm), we obtain 



_d_ 

dm 





d'tt 2 Q4tt 4 A d' 
F + 



dm m 2 M 2 dm \A 



G\ 64vr 4 A d 



2„2\ 



M 2 dm 



v 



0. 



(24) 



with the derivative d' /dm acting only on the integrals of motion. Finally, we use "Bernoulli's 
relativistic equation" 7 2 — u 2 = 1, which, given the definitions of the integrals of motion E(ty) 
and L(ty), can be written as 

/ Art 1 \ 

A 2 




e 2 x 2 (A-M 2 ) x 2 E 2 n 2 (0)L 2 x 2 A 2 



fi 2 T] 2 



M 4 



+ 



fl 2 7] 2 



fI 2 T] 2 



M 4 



where we introduced the dimensionless variables 



(25) 



x 

y 



tt F (0)m 
a^/^Q. 



Vim , 



As a result, substituting the right-hand part of (^) into the first term of 
differentiation, we obtain the first first-order differential equation 



H 2 i] 2 



+ 



n 2 (o) 



x 2 - 1 



M 2 dy 
~dx 



dM 2 M 6 Ql{0)L 2 xM 2 



(26) 
(27) 

and performing 



dx 
1 



x 3 A 



jj?7] 2 



A Q 2 (0) \fi 2 r] 2 



-2A 



de 2 



+ 



x 



dtt 2 



fi 2 V 2 dy Q 2 {0) dy 



-2 1 
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The second first-order differential equation is Bernoulli's equation (|25|), which should now be 
considered as an equation for the derivative dy/dx. The set of equations (^) and (p8|) allows a 
general solution to be constructed for a one- dimensional jet immersed in an external magnetic 
field. 

We emphasize one important advantage of the set of first-order equations ( ^5|) and ( |28|) over 
the initial second-order equation fllQl). The point is that the relativistic equation (0), which is 
basically the force balance equation, contains the electromagnetic force 

F em = p e E + jxB, (29) 

in which the electric and magnetic contributions virtually cancel each other far out from the 



rotation axis w 3> Rl- Using Bernoulli's equation (^), we can derive [[31 

|p c E + jxB| 1 



|j x B| 7 



2 ' 



(30) 



When analyzing (|T0"D, we therefore must retain all higher order terms ~ j , while the zero- 
order quantities p e E and j x B in (|2~8D are analytically removed using Bernoulli's equation, so all 
terms of this equation are of the same order. 

Finally, it is also important that the exact equation (|2~8D has no singularity near the rotation 
axis. In other words, its solution contains no 5-shaped current / oc 5{w) flowing along the jet 
axis; several authors pointed out to the necessity of it |7|, |32fl . 



3 Exact solutions and numerical results 



Let us consider the basic properties of the set of equations ( |25| ) and (^). As can be easily verified, 
in the relativistic case under consideration, we may assume 7 = u z with high accuracy. Far from 
the rotation axis, x ^> 7i n [w ^> 7 m -RL), equation ( p5|) can be rewritten in the limit M 2 -C x 2 as 

dV _ 8n 2 E(^) 

dzu zun 2 (^) { ' 

or, equivalently, 

As we see, equation (|3l|) does not contain M 2 at all and can therefore be integrated indepen- 
dently. This must be the case, because equation fl3~T|) must coincide with the asymptotics of the 
force-free equation, which can be derived from (|25|) by going to the limit M 2 — > 0. Assuming now 
that B z (rj) = B cxt in (|32|) , we obtain, in particular, for the jet radius 



r 2 - lim AnEW (33) 



Consequently, the jet radius is determined by the limit of the E(^f)/Qp(^f) ratio as ^ — ► 
In particular, for Efif) and fip^) given by (|18|)- ([20|) , we have 



E(^) 1 

lim _\ \ = —^0, (34) 



7 



so the limit (^) does actually exist. As a result, we obtain 



n = i/=£-. (35) 



7T 



which essentially coincides with estimate (|2]). This is no surprise, because we show below that 
equations ( p5|) and (p8|) for the integrals of motion (|18D, (p0|), and ( p3|) have a constant magnetic 
field as their solution over a wide range of w. 

Let us now consider in more details the behavior of the solution in the inner jet region, where 
\& \l/o an d, hence, the integrals of motion can be approximately written as 

m = (36) 

f2 F (tf) = tt = const, (37) 
r](^) = T] = const, (38) 

with E(ty) = 7i n /i77o + ^f-^ and e(\I/) = 7^770 = const. As a result, QpL/ [irjo = 2y, and we can 
rewrite equations (|2"5|) and (EBI) as 



J _ m 2 ) 2 ( ± V = %Hl - x 2 - 2M 2 ) + a; 2 ( 7in + 2y) 2 - % 2 - ^(1 - x 2 - M 2 ) 2 , (39) 



/ 2 2 J^ 2 „ ,,2 7in^ 2 4|/ 2 M 6 

^ + x - »-*r = 2lM - (i - ^ - m 2 ) + xH i =p - a/ 2 ) • < 40 » 

Equations (^) and (|40f) describing the internal jet structure can be solved analytically. It can 
be verified by direct substitution that we have the following asymptotics for x <C 7i n : 

M 2 (x) = M 2 = const, (41) 

V(x) = (42) 
which correspond to a constant magnetic field 



47r7 in /i^ 7i 



B z = BM = = ^B(R L ) = const, (43) 

where B(R^) = ^o/Rl- Here, we assume that 7; n ^> 1, which is the typical for jets from AGNs 
and radio pulsars. 

The solution of fl39"D and fllPf) for x ^> 7 in , i.e., at w jinRh, depends on the relationship 
between 7 in and M = M(0). For example, at Mq > 7 2 n , when, according to ([|3|), the axial 
magnetic field is fairly weak, the total magnetic flux within w < 7i n i?L 

^(7ini? L ) « 7r&R£fl,(0) (44) 

can be written as 

^(7ini?L) ~ ||*in, (45) 
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where the flux \l/ m is given by (^). We see that, if the condition M5 2 > 7 m is satisfied, then the 
total magnetic flux within w < 7i n -RL is lower than \l/ in ; so, outside this region, the particles also 
make the main contribution to E(ty) as before, while the contribution of the electromagnetic field 
may be neglected. As a result, at w ^> 7i n -RL> the solution of (^) and ( [40"D has a quadratic rise 
of M 2 and a power- law falloff of the magnetic field [^, ^3|, |24| : 

M 2 {x) = M 2 ^ > x 2 , (46) 

Tin 

B s (x) = B g (0)^. (47) 

Consequently, the magnetic flux increases very slowly (logarithmically) with the distance from 
the rotation axis: 

^(x) oc ln(x/7 in ). (48) 

Such a behavior of the magnetic field, in turn, shows that the transition flux \I> = \l/ m is reached 
exponentially far from the rotation axis, in conflict with the estimate (§) corresponding to the 
assumption of jet collimation. 

We may thus conclude that an external constant magnetic field limits the Mach number at the 
rotation axis above 

M 2 (0)<M 2 ax = 7m . (49) 



Accordingly, as follows from (f43|), the magnetic field at the rotation axis cannot be weaker 
than 

£Un = —B(R L ). (50) 
0"7in 

If, however, the Mach number at the rotation axis does not exceed 7; n (i.e. if Mq < 7 m ), then, 
as for the similar asymptotics w <^ ji n RL, the solution of (|39"D and (|40[) for 7i n -RL ^ w *C rj gives 
a constant magnetic field (£H), which corresponds to the solution 

y(x) = 2M| X2 - (51) 
At the same time, in this case, we have only a linear increase in the square of the Mach number 

M 2 (x) = M 2 — < x 2 . (52) 

7in 

Then, according to (p7|) and (|51"D, the jet radius can be written as 



r 



" MS R L , (53) 



which is equivalent to fl35|) [and in agreement with (Q)]. Moreover, as can be easily verified, the 
constant magnetic field B z = B(0) for the invariants (p~8|) (pQ|) proves to be an exact solution of (|3l| ) 
in the entire jet up to the jet boundary, w = r- y Here, we may therefore assume -B(O) = -B ex t- 
Consequently, according to ([43D , we obtain 

Ml = *™ (54) 

O" -Dcxt 
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Using relation (|54|), we can also express all the remaining jet parameters in terms of the external 
magnetic field. 

Note that the absence of a declining solution B z oc zu~ 2 [see (|)] is associated with the first 
term in the right-hand part of ( p8|) proportional to L 2 . This term, which changes appreciably the 
behavior of the solution, appears to be missed previously. As it was already emphasized above, 
this is not surprising because the corresponding term in the second-order equation (pi]) is of high 
order and small. On the other hand, far from the rotation axis w ^> j- m , equation (^) can be 
rewritten as 

d 



dx 



M 2 \ VL Fm x%x 2 + M 2 ) ' 1 ' 

in which both terms are of the same order. Neglecting the term proportional to L 2 , we arrive at 
the solution (fl6|), M 2 oc x 2 , for Qp = const and 77 = const. The conservation of function 

HT]Q F x 2 

H = -JP- (56) 

was first found by Heyvaerts and Norman f?J for conical solutions, when all quantities depend 
only on spherical coordinate 9, but has also been repeatedly discussed when analyzing cylindrical 
flows. However, as we see, H is generally not conserved in the cylindrical geometry for relativistic 
jets. To be more precise, the second term in (|55|) turns out to be significant for all models with 
a nearly constant density of the longitudinal electric current in the central jet region, where the 
invariant L(ty) linearly increases with magnetic flux \1/ if \1/ -C ^o- 

Thus, we conclude that the solution with a central core (Q) cannot be realized in the presence 
of an external medium with a finite regular magnetic field. This conclusion appears to be also 
valid in the presence of a medium with finite pressure P. Indeed, since the magnetic flux (]48| ) 
increases very slowly (logarithmically), the solution @ yields an exponentially large jet radius 
rj ~ i?L exp (^o/^m)- Accordingly, the magnetic energy density must also be low at w ~ rj. 
However, this configuration cannot exist in the presence of an external medium with finite pressure 
P, irrespective of whether it is produced by a magnetic field or by a plasma. We may therefore 
conclude that the solutions with a central core can be realized only for a special choice of the 
integral L(^), which increases only slightly with the magnetic flux, and only in the absence of an 
external medium. For the most natural (from our point of view) models with a constant current 
density in the central jet regions, the solution with a central core cannot be realized even in the 
absence of an external medium. 

In order to derive now the energy distribution in the jet and the particle Lorentz factor, it is 
convenient to introduce the quantity 

g(x) = — (57) 
x z 

Since at large distances w ^> 7i n -RL, according to (0), we have 

7ffl 9 



E g+V 



(5f 



g is simply the ratio of the energy flux transferred by particles W pai t to the energy flux of the 
electromagnetic field. As a result, from relation fl52l ) for x ^> 7^, we obtain 



^ - ^ « 1. (59) 

Wtot 7in 
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Accordingly, from (0) and fl5"8|) for the particle Lorentz factor at x ^> 7^, we derive 

j(x) = x. (60) 

Finally, expression ( |l3|) for the (^-component of the 4-velocity vector u v yields the following 
toroidal velocity v v = u ip / r y at x ^> 7 in : 

^(x) = -. (61) 
x 

We see that the particle energy approaches the universal asymptotic limit (|(J) at w ^> 7i n -RL- 
Naturally, such a simple asymptotics can also be derived from simpler considerations. Indeed, 
using the "frozen- in" equation E + vxB = 0, we obtain for the drift velocity 

u «=w = (w + w) 1 - (62) 

In our case, however, according to (ffl) and ([?]), we have 

B% » |E| 2 , (63) 

|E| 2 ^x 2 5 2 . (64) 

As a result, relations d62|) -(|64D immediately lead to the exact asymptotics (|60|) . It thus follows 
that, for example, for electron-positron jets from AGNs, the jet particle energy is typically 

E « m c c 2 ^- ~ (10 4 - 10 5 ) Mev. (65) 
Rl 

On the other hand, according to ([57]), we reach a very important conclusion that, far from the 
light cylinder, w 7i n -RL, 

g < 1. (66) 

Consequently, according to (|58|), the particle contribution to the general energy flux balance 
proves to be minor. For example, at B cxt ~ -B m in for r ~ rj, we have 



Wpart /77 



While in the general case, we obtain 

~B{R h ) 



(67) 



Wpart 1_ 
W tot ~ a 



1/2 

(6£ 



Thus, we reach a fairly nontrivial conclusion that the fraction of energy transferred by particles 
in a one- dimensional jet must be determined by the parameters of the external medium. 

Let us now discuss the results of exact calculations obtained by numerical integration of equa- 
tions ( j25| ) and ( |28| ) with the integrals of motion (]T8|) (pQ|) and (^3|). In figures [L](a) and [j](b), the 
Mach number and the energy flux concentrated in particles 7/^7 are plotted against x = VL-^vj for 
Mq = 16, 7i n = 8, and a = 1000. The dashed lines indicate the behavior of these quantities that 
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follows from the analytical asymptotics ([52]) and fl60|). As we see, at sufficiently small x when the 
integrals of motion fll8|)-(|2(^) and ([23]) are similar to (|36|)-(|38D, the analytical asymptotics match 



the exact numerical results. On the other hand, as expected, 7/^77 and M are zero at \I/ = \&0) 
i.e., at the jet edge. Figure |l|(c) shows the dependence of the poloidal field, x~ 1 (dy/dx), for the 
inner parts of the jet, $ < ^o- We see from this figure that the magnetic field is nearly constant 
at x > 7; n , in agreement with the analytic estimates fl43|) and flST]). Of course, in general, the 
structure of the poloidal magnetic field is determined by a specific choice of the integrals E(^f) 
and 

In conclusion, it is of interest to compare the energy of jet particle with the limiting energy 
acquired by the particles as they outflow from the magnetosphere of a compact object with a 
monopole magnetic field. According to calculations by Beskin et al. ||33|| , the particle Lorentz 



factor outside the fast magnetosonic surface r > <r 1//3 i?L in the absence of an external medium can 
be written as 

7(1/) = y 1/3 , y>7in, (69) 
7(2/) = Tin, V<lL (70) 

where y is given by (^). On the other hand, relations ( |HT| ) and (^) for the jet yield 



*>-(*)>. y>p 

i{y) = 7m, y<-jj|. (72) 

As shown in figure ||(a), for M 2 > 1, i.e., for B cxt < B CT , where 

Bcr = ^B(R L ), (73) 
cr 

the Lorentz factor of the jet particles ([H]) is always larger than the Lorentz factor acquired by the 
particles as they outflow from a magnetosphere with a monopole magnetic field, but, of course, is 
always smaller than the critical Lorentz factor 

i(y) = y, (74) 

which corresponds to the complete transformation of the electromagnetic energy into the particle 
energy. This implies that, at B cxt < B CT , the particles must be additionally accelerated during the 
collimation coupled with the interaction of the outflowing plasma with the external medium. If, 
alternatively, B cxt > B CT , then, in the inner jet regions, at 

- < ||^. (75) 

the particle energy on a given field line turns out to be even lower than that for a monopole 
magnetic field, as shown in figure 0(b). 

The latter result can be easily explained. Indeed, for the integrals of motion (|3T)|)-(|3"BD we 



consider, the factor D, whose zero value determines the location of the fast MHD surface (see [27| 
for more detail), can be rewritten in the case of a cold plasma as 

B 2 Ax 2 v 2 M 4 

^ d = a+ 4= a + wm^w^ - < 76) 
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Figure 1: The dependences of Mach number M 2 (a), energy flux 7/xr/ concentrated in particles 
(b), and poloidal magnetic field B z (c) obtained by numerical integration of the equilibrium 
equations ( p5|) and (p8|) for M = 16, 7i n = 8, and a = 1000. The dashed lines indicate the 
behavior of these quantities that follows from the analytic asymptotics. 
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ig(y) 




Figure 2: Logarithmic dependence of the particle Lorentz factor on the magnetic flux for a 
monopole field 7 ~ y 1 ^ 3 , for a one- dimensional jet 7 ~ (MQ/7; n ) 1 / 2 ?/ 1 / 2 , and for complete transfor- 
mation of the electromagnetic energy into the particle energy 7 ~ y at: (a) M 2 > 1 (B ext < -B CT ); 
and (b) M 2 < 1 (-B ex t > B CT ). The intersection of the straight lines (M 2 pa 1) determines the 
position of the fast MHD point. 
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It is easy to show that expression ( |76|) for y and M 2 given by ([51]) and ([52]) is negative at 



M 



(77) 



i.e., at x corresponding to ([T5|). Consequently, we may reach another important conclusion that, 
for sufficiently strong external magnetic fields B ext > B CT (^) when M 2 < 1, a region with a 
subsonic flow inevitably emerges in the inner jet regions w < r s , where 



a 



B, 



ext. 



B(R L ) 



3/2 



ft 



(7c 



At the same time, a region with the subsonic flow can be produced far from the compact object 
either by a shock wave or by a strong distortion of the magnetic field within the fast magnetosonic 
surface located in the vicinity of the compact object. In both cases, the magnetic-field perturbation 
causes the particle energy to decrease. 



4 Discussion 

Thus, we conclude that the exact equilibrium equations fl2"5|) and (|28| ) do actually allow a con- 
struction of a self-consistent model for a jet immersed in an external uniform magnetic field. The 
advantage of these equations over equation fllQ|) results from the fact that all terms in ( p8|) are 
of the same order. In this case, the uniformity of the poloidal magnetic field within the jet (|43| ) 
results from the choice of integrals (|5o|)-(|5£|). In general, the poloidal magnetic field depends on 
the specific form of the integrals. A full analysis of the possible solutions is beyond the scope of 
this paper. 

We have shown that the fraction of energy transferred by particles W part /Wtot must be largely 
determined by the parameters of the external medium. In case o » a CI , where 

B '^l 1/2 , (79) 



B 



ext. 



the energy transferred by particles is only a small fraction of the energy flux W em transferred by 
the electromagnetic field. Consequently, the jet is strongly magnetized (W^art "C Wtot) only at 
sufficiently large a. If, however, the magnetization parameter does not exceed <r cr , then, in this 
case, an appreciable part of the energy in the jet is transferred by particles. This, in turn, implies 
that a considerable part of the energy must be transferred from the electromagnetic field to the 
plasma particles during jet collimation. It is interesting that cr cr turns out to be approximately 
the same both for AGNs and for fast radio pulsars: 

cr cr « 10 5 - 10 6 . (80) 

We have shown that the central part of the jet must be subsonic for sufficiently strong ex- 
ternal magnetic fields. Thus, the theory gives direct predictions whose validity can be verified 
by observations. It should also be noted that the results obtained above are applicable both to 
electron-positron and to electron-proton jets. However, a direct evidence that the jets in AGNs 
are actually electron-positron ones has recently appeared [|33], ft5j] . 
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In our view, an important result is that, if the external regular magnetic field is taken into 
account, the MHD equations allow a self-consistent model to be constructed for a jet with a 
zero total longitudinal electric current, /(^o) = 0. In this case, a uniform magnetic field that 
matches the external magnetic field can also be a solution for the inner jet regions. As was 
already emphasized above, the radii of the jets from AGNs can thus be also explained in a natural 
way. In addition, since only a small fraction of the electromagnetic-field energy is transformed 
into the particle energy, the energy transfer from the compact object in the region of energy 
release can be explained as well. At the same time, extending the MHD solution to the jet region 
requires very high particle energies (~ 10 4 MeV), which have not been recorded yet. However, a 
consistent discussion of the outflowing-plasma energy requires a proper allowance for the particle 
interaction with the surrounding medium (for example, with background radiation), which may 
cause a significant change in particle energy. 

As for the quantitative predictions about the real physical parameters of jets, they, as we 
showed above, essentially depend only on the following three quantities: the magnetization pa- 
rameter cr (|17D, the Lorentz factor 7; n in the generation region, and the external magnetic field 
B ext . In this case, the main uncertainty for electron-positron jets from AGNs (Bi n ~ 10 4 G, 
i?i n ~ 10 14 cm) is the value of the magnetization parameter. Indeed, this quantity depends on 
the efficiency of pair production in the magnetosphere of a black hole, which, in turn, is deter- 
mined by the density of hard gamma-ray photons. As a result, if the density of hard gamma-ray 
photons with energies J5 7 > 1 MeV near the black hole is high enough, then the particles will be 



produced by direct collisions of photons 7 + 7 — > e + + e~ . This causes an abrupt increase 
in the multiplicity parameter A = n/ncj ~ 10 10 -10 12 , where ncL ~ QB /2nc is the characteristic 
particle density required to shield the longitudinal electric field. Using the well-known estimate 
(see, e.g., 

neB m Rf n 

a ~ — — , (81) 

Xm e c 6 

we obtain 

a ~ 10 - 10 3 . 7in ~ 3 - 10. (82) 

On the other hand, for low densities of gamma-ray photons when an electron-positron plasma 
can be produced only in regions with a nonzero longitudinal electric field, which are equivalent 
to the outer gaps in the magnetospheres of radio pulsars [[T7|, |58fl , the multiplicity of the particle 
production is fairly small: A ~ 10-100. In this case, we obtain 

a ~ 10 11 - 10 13 , 7 in ~ 10. (83) 

Finally, for fast Crab- or Vela-type radio pulsars {B m ~ 10 13 G, polar-cap radius R\ n ~ 10 5 cm, 
and A = n/riQj ~ 10 4 ) in which jets are observed, we have [31] 

a ~ 10 6 - 10 7 , 7 in ~ 10 2 - 10 3 . (84) 

Relations (jS^ ) and (jB3|) show that the properties of the jets from AGNs considerably depend 
on the magnetization parameter cr. For example, according to (|79|) , the jet particles for sources 
with large a ~ 10 12 transfer only a small fraction of the energy compared to the electromagnetic 
flux, so the flow within the jet differs only slightly from the force-free flow. In addition, in this 
case the external magnetic field B ext ~ 10 -6 G exceeds the critical magnetic field B CT ~ 10~ 7 G. 
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According to ( f43| ) and (|73|), this implies that a subsonic region must exist in the inner regions 
of such jets. On the other hand, a substantial part of the energy in sources with magnetization 
parameter a ~ 100 during jet collimation must be transferred by plasma particles, and no subsonic 
region is formed near the rotation axis. As for the fast radio pulsars, the condition 5 ext <C B CT is 
satisfied for them, so a subsonic region in the central parts of pulsar jets is not achieved either. 
On the other hand, the estimate (|79|) shows that an appreciable part of the jet total energy must 
be coupled with particles. 



We emphasize that, since the jet radius fl53|) for AGNs always exceeds the light-cylinder radius 
by several orders of magnitude, 

rj ~ (10 4 - lO 5 )^, (85) 

the toroidal magnetic field within the jet must exceed the poloidal magnetic field B p in the 
same proportion, 

B v « ^-B p ~ (10 4 - 10 5 )S cxt . (86) 

Consequently, detection of such a strong toroidal component would be a crucial argument for 
the picture discussed here. Unfortunately, determination of the actual physical conditions in jets 
currently involves considerable difficulties. Nevertheless, not only data on the direct detection 
of such a structure |39[] but also evidence for the existence of magnetic fields B ~ 0.1 G, closely 



matching the estimate (|86|) [p4| , have recently appeared. 
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